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NOMENCLATURE 


f = local fraction of inner jet fluid (ie., concentration) 
k = turbulence energy [m^/s^] 
m = mass flow rate [kg/s] 
r = radius (mm) 

Rq = radius of sudden expansion (mm) 
u = axial velocity (m/s) 
v = radial velocity (m/s) 
w = azimuthal velocity (m/s) 
z = streamwise coordinate (mm) 

= dissipation rate of turbulence [m^/s^] 
m = eddy diffusivity 


Superscripts 

mean quantity 
' = fluctuating quantity 
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INTRODUCTION 


The outlet flow for a combustion chamber is determined by the 
complex pattern of the flow field through it. In order to perform 
parametric studies on the flow, it is necessary to have an 
accurate simulation program for prediction of the combustion 
process. Since combustors of practical interest are highly 
turbulent, understanding turbulent transport in this type of flow 
is critical to the development of computational procedures to be 
used in these simulations. 

Modern combustion chambers typically have an annular 
configuration with fuel in a fine spray form being mixed into the 
airstream. This mixture is achieved with an array of fuel nozzles 
which ensure uniform fuel spray distribution. these nozzles may 
be of the swirl or non-swirl type, although the swirl type is 
usually used to promote rapid burning. The fuel-air mixture is 
burned in the combustion chamber and the resulting gases are 
delivered to the turbine inlet. A stabilized flame in the 
combustion process is necessary to achieve uniform burning, 
resulting in a desirable temperature profile at the turbine inlet 
(pattern factor) . Other important parameters in the combustion 
process include burning length (flameout), noise production, 
pollutant emissions (reactant products), and engine performance 
(combustion efficiency). 

Computational procedures for predicting the combustion process 
are being developed and improved by numerous researchers (1,2,3,). 
These procedures predict the velocity, species, temperature, and 
reaction rate distribution within the combustors which in turn are 


used to calculate the previously named parameters. Because of the 
turbulent nature of the flow, mathematical models are used for the 
turbulent transport of mass (species), momentum, and heat. The 
data used to develop and verify these models have in the past been 
restricted to velocity and momentum transport measurements. 

Methods used prior to the mid-nineteen seventies for acquiring 
turbulent mass and momentum transport data have been indirect, 
requiring compromising assumptions or probes unsuitable for 
recirculating flows. Bennett & Johnson (4) have used new optical 
techniques to simultaneously measure a scalar quantity (mass) and 
velocity, and therefore mass transport. 

Morganthaler (5) performed an 

"assessment of the relative importance of turbulent 
mass, momentum, and energy transport, so that emphasis 
for the future could be directed to modeling the cri- 
tical processes rather than merely continuing the 
historical trend of modelling turbulent momentum 
transport . " 

His results point out a critical need for the future in that 

"turbulent transport of mass was demonstrated to be far 
more significant than the transport of either energy or 
momentum for a coaxial hydrogen jet reacting with an 
external high temperature air stream." 

In order to provide an adequate data base for verifying 
computational procedures used for modeling mass transport, Bennett 
& Johnson (6) used laser velocimetry (LV) to obtain velocity 
measurements while simultaneously using laser-induced flourescence 
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(LIF) to obtain concentration data for a coaxial jet geometry. In 
a companion paper, Syed & Sturgess (7) used the data of Bennett & 
Johnson for comparison with computational methods currently in 
use. While finding qualitative agreement of all variables, there 
were some significant quantitative discrepancies. Among their 
conclusions was that "prediction of mean quantities alone is not a 
sufficiently strict criterion" for evaluating the success of a 
turbulence model, but that "it is necessary to examine also the 
agreement obtained for fluctuating quantities and their cross 
correlations" (i.e., mass flux). 

It should be pointed out that some of the difficulty 
encountered in accurately predicting the flow can be attributed to 
the lack of a set of good initial conditions. Another limiting 
factor - possibly the dominant one in certain situations - is the 
incomplete understanding of the mechanisms involved in the flow. 
One such mechanism is the influence of the large-scale motion on 
the turbulent transport. The understanding of this influence is 
the main thrust of this research effort. 

Other comparison studies of turbulent flow in a coaxial jet 
have been done by Habib & Whitelaw, including experimental data 
obtained using hot wires (8) and using LV (9). In each study, the 
largest discrepencies were in the upstream region where the flow 
is developing and recirculation zones occur. They concluded that 
these errors were "associated with the inadequacy of the 
eddy-viscosity hypothesis." Again it is seen that a better 
understanding of the mechanisms involved is needed to improve the 
computational procedures currently in use. 


The experimental efforts of Bennett & Johnson have also 
continued and produced some interesting results (10) which in some 
cases are contradictory to traditional thinking. Complete details 
of the experiments and data are found in the final NASA contractor 
reports of Bennett & Johnson (11) for a non-swirling flow and 
Roback & Johnson (12) for a swirling flow. These works were done 
on a coaxial jet configuration using water in both jets and 
fluorescin dye as a trace element. Bennett & Johnson show two 
distinct shear regions for the non-swirling flow (fig. la); one 
between the inner and annular jets and one between the annular jet 
and the recirculation region. The innermost shear layer develops 
as the annular fluid gradually fills the center jet, resulting in 
what has been called counter-gradient transport. Computational 

procedures for evaluating scalar (concentration) gradient and a 

£ 

transport diffusion coefficient, ni defined by 

m - £ m /1I\ 

In many turbulent flows the scalar transport does not follow such 
a simple dimensional model and the traditional approach does not 
work . 

The swirling flow results of Roback & Johnson (12) show a flow 
field (Fig. lb) similar to that for the non-swirling case, they 
found a large eddy shear region between the inner and annular jets 
and one between the annular jet and outer recirculation zone. 
Unlike the non-swirling flow however, they found a large 
recirculation region on the centerline which gave the flow some 
different characteristics. 
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The multiple scale results of Bennett & Johnson are an example 
of a situation for which this simple gradient model is inadequate. 
Figure 2 shows regions of the flow field where this is the case. 
This region is qualitatively similar to the region where the inner 
jet fluid is being accelerated by the annular jet, although it is 
somewhat more extensive. Bennett & Johnson hypothesized that this 
was due to "response time or distance required to change the 
character of the turbulent structure" . Another interesting 
finding was that the peak absolute values of the axial mass 
transport rates were higher than the values for the radial 
transport even though the peak radial concentration gradients were 
more than ten times the axial concentration gradients. 

Observation of the concentration signal during the experiments 
showed much large-scale flow intermittency with the developing 
portion of the flow field. At axial locations within the upstream 
region (where the flow is developing), there were large-scale 
fluctuations in the signal: they were either negative (indicating 

a slug of annular jet fluid) or positive (suggesting a slug of 
inner jet fluid) . No immediately recognizable periodicity in the 
occurrence of these fluctuations was found and no local peaks in 
the concentration autocorrelations were identified. It was 
believed that these excursions did identify the presence of a 
large-scale motion within the region. 

This later work of Bennett & Johnson gives some justification 
for studying the influence of the large-scale motion on the 
turbulent transport. The comparisons of Syed & Sturgess, as well 
as Habib & Whitelaw, also support this in that, for both studies. 


the region of greatest discrepancy is the upstream region where 
the flow is developing. Since this region is known to include 
large-scale motion, it seems reasonable to investigate the effects 
of that motion on the flow. 

In addition, Schetz (13) has reviewed the mixing flows in 
terms of both experimental results and the implications for 
turbulence modeling. One of his basic conclusions was that 
large-scale structures are an essential aspect of future 
experimental and modeling studies for a wide variety of flows. 
Mathiew & Jeandel (14) point out that the large interacting eddies 
can have a time scale quite differenct from the dissipative time 
scale associated with the fine structure portion of the flow. 
Launder (15) has introduced multiple time scale models, chosen to 
give the best predictions. These scales could be chosen more 
appropriately from experimental data directly if the role of the 
large scale were identified. Finally Borghi (16) questions the 
importance of large coherent structures in combusting flow. His 
conclusion is that, for the upstream region, the numerical 
predictions are not close to observed data because either large 
structures or the redistribution of kinetic energy is badly 
calculated . 

In summary, for the confined coaxial jet, the presence of 
large-scale structures has been identified in the region where 
computational predictions have shown to be least effective. 

Several researchers have hypothesized that the large-scale 
structures may be the source of these discrepancies. With this in 
mind, this study was undertaken to identify a way of detecting 
these structures and to determine the influence they have on the 
turbulent transport of the flow. 


BACKGROUND 


The notion of large-scale coherent structures in turbulent 
flow is not a new one. Several authors - including Roshko (17), 
Cantwell (18), and Davies & Yule (19) - have summarized many of 
the developments in this area. A starting point in any discussion 
must be a clear definition of the subject. Yule (20) defines 
coherent structures as large eddies which, (i) are repetitive in 
structure, (ii) remain coherent for distances downstream very much 
greater than their length scales, and (iii) contribute greatly to 
the properties of turbulence, in particular, turbulent energy and 
shear stress, entrainment and mixing. 

The experimental work of Brown & Roshko (21) is considered 
classic in the study of large-scale motion in turbulent mixing 
layers. Their efforts (fig. 3) revealed the "presence of 
well-defined large structures", with the Reynolds number varying 
from a low value with no visible fine-scale turbulence to a higher 
value where it does exist. The significant point is that in each 
case "the measured mean properties of the flow, the velocity and 
density profiles, spreading rate, etc. are the same". They 
concluded that the "mean flow is controlled by the large organized 
structures which, it may be seen, are not affected by the 
small-scale turbulence appearing at the higher values of Reynolds 
number . " 

By observing the movement of the structures. Brown & Roshko 
(21) found that their spacing and diameter increased with 
increasing downstream distance. The large-scale structures move 


at nearly constant streamwise velocity (equal to the mean flow 
velocity) which is independent of their size and location. Brown 
& Roshko conclude that they amalgamate into larger structures as 
they convect downstream. Winant & Browand (22) observed fluid 
rolling up into "discrete two-dimensional vortical structures" 
which periodically "interact by rolling around each other" 
eventually forming one larger vortex. Finally, Dimotakis & Brown 
(23) found similar large-scale coherent structures at a high 
Reynolds number (3x10^). Their flow visualization study showed 
that these structures did exist and "appear to dominate in 
determining the overall characteristics of such flows". 

Blackwelder & Kaplan (24) studied bursting in a turbulent 
boundary layer. They observed, using flow visualization, "a high 
degree of coherence over a considerable area in the direction 
normal to the wall", and found that bursts were associated with a 
high degree of velocity fluctuation. Their conditional averaging 
process showed the coherent structures to have turbulent transport 
properties an order of magnitude greater than the overall flow. 
This is in agreement with the work of Lu & Willmarth (25) who 
found large contributions to the turbulent transport from the 
coherent structures. Boundary layer .flows have been studied more 
than any other type and are probably the best understood in terms 
of coherent structures. It is useful to apply this knowledge in 
developing a better understanding of other flows. 

The experimental work of Bennett & Johnson (11) for 
non-swirling flow and Roback & Johnson (12) for swirling flow is 
the basis for the present study. The configuration they chose was 

a confined coaxial jet with a sudden expansion, a meaningful yet 

8 


geometrically simple case. It simulates the two-stream inlet and 
the diffuser (sudden expansion) similar to flow in combustors. 
Water was chosen for the experiments; the resultant capability for 
injecting dye into either stream was utilized for flow 
visualization and mass transport studies. The use of water limits 
the simulation in that it eliminates combustion and multiple 
species as variables to be considered. Water also reduces the 
effect of molecular diffusion to a negligible level. These 
limitations still resulted in a set of results useful in 
evaluating combustor design calculation codes; before the addition 
of combustion, the codes should be expected to accurately predict 
this more limited flow. To create this data base for comparison 
with models, the emphasis was on statistically repeatable data 
collection. As the emphasis was not on large-scale structures 
(i.e., conditional sampling), this limits the data's usefullness 
in studying large-scale structures somewhat. For example, it is 
not possible to trace the flow of a particular structure past a 
measurement location since the data collection rates needed for 
statistically steady data limit the number of samples within any 
one structure. 

Bennett & Johnson (4, 11) and Roback & Johnson (12) documented 
the flow fields shown in figure la and lb respectively. The 
ratios of jet diameters are approximately 0.25 for the inner jet 
and 0.50 for the annular jet (compared to the pipe diameter). The 
flow chosen for documentation has Reynolds numbers of 15,900 and 
47,500 for the inner and annular streams respectively. These 
values are within the turbulent flow range and are typical of gas 
turbine combustors. The data base resulting from these 
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experiments consists of simultaneous two-component velocity data 
(u-v, u-w, v-w) and simultaneous concentration-velocity data (u-f, 
v-f, w-f) at numerous axial and radial locations on the flow 
field. 

The comparison work of Syed & Sturgess (7), which utilized 
Bennett & Johnson's data base, was done with a two-dimensional k-e 
model employed in the TEACH computer program. The largest of the 
discrepancies between experiment and predictions occurred in the 
upstream region where the flow was developing and included 
recirculation zones. For fluctuating velocity components, 
agreement was good qualitatively, but only fair quantitatively. 
Mean axial velocity profiles were found to be in better agreement. 
The centerline concentration growth is predicted to occur sooner 
than experiments indicated, while the subsequent decay is 
predicted later than found experimentally. Mass flux predictions 
exhibited similar discrepancies, as did limited momentum transport 
comparisons . 

In the two Habib & Whitelaw studies (8, 9), their data is 
compared with predictions, again using the two-equation turbulence 
model in the TEACH code. Using hot-wire data, predictions of mean 
velocity distributions agreed well with the experimental data in 
the downstream regions. Shear stress distributions showed good 
agreement downstream but again had some trouble in the upstream 
flow developing region. Although the LV measurements in the 
second study were an improvement over the hot-wire data, errors 
were still found in the recirculation zone, reinforcing the need 
for modeling improvements in the prediction of turbulent 
transport . 
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In summary then, large-scale structures have been shown to 
contribute significantly to turbulent transport (at least in 
boundary layers). In addition, flow regions for which numerical 
predictions and experimental results differ the most are known to 
contain large-scale structures. It seems apparent that a better 
understanding of large-scale structures will be necessary for 
continued development of prediction methods. The present 
investigation represents an initial effort in this area for 
axisymmetric shear layers. 
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ANALYSIS 


This study involved looking at the effects of the large-scale 
structures on the data collected by Johnson St Bennett (11) and 
Roback & Johnson (12). These data were used previously to study 
the overall flow characteristics. The present effort verified the 
overall properties; but the main thrust was on the large-scale 
structure properties and their influence on the overall flow. 

In order to identify the role of large-scale structures in 
influencing turbulent transport, it is first necessary to separate 
these structures from the small-scale turbulent portion of the 
flow. A detector or trigger for the structure must be selected. 
Several potential detectors were identified based on other 
studies, including several done on turbulent boundary layers (24, 
26). One possible detector is large velocity fluctuation compared 
to turbulence intensity, the rationale being that large-scale 
structures will involve blobs with larger velocity changes than 
for the overall flow. Another possible detector is concentration 
fluctuation. The entrainment and vortex pairing phenomena 
previously discussed led to this idea. Roshko (17) showed data 
from a concentration probe in a mixing layer where fluid from one 
stream penetrates deeply into the other stream, resulting in large 
concentration fluctuations. He states that information about the 
mixing can be obtained by measuring a scalar property such as 
temperature, density, or species concentration. The experiments 
of Bennett & Johnson and Roback & Johnson, previously discussed, 
also showed these large concentration fluctuations in regions 
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where the large-scale structures are known to exist. 

Concentration fluctuation was chosen as the detector for this 
study . 

The ultimate test of any detector is the completeness of the 
emerging pattern. Completeness refers in this case to a 
well-defined pattern, since a detector which yields randomly 
varying results cannot be a true detector. One way of testing a 
detector is to plot the relative frequency of each individual 
detection as a function of position. Relative frequency is the 
number of occurrences at each location ratioed by the complete set 
of data at that location. This can help identify the approximate 
size and most probable location of the large-scale region. These 
plots were done and will be discussed later. 

Another useful plot is the scatter plot where the velocity 
fluctuation is plotted as a function of the concentration 
fluctuation or where fluctuation of one velocity component is 
plotted against the fluctuation of another velocity component. 
These types of plots help establish the existence of the large 
scale in the flow and assist in separating these structures from 
the general low level turbulence. They are also beneficial in 
determining if the large-scale correlation is significant. 

An example of a plot with a zero correlation, shown in Figure 

4, has almost all of the data concentrated uniformly near the 

center. A small amount of data has fairly large negative 

concentration fluctuations, but are equally distributed among 

positive and negative velocity fluctuations, yielding a zero 

correlation. Figure 5 shows data with much more scatter in both 

directions and a zero correlation. A situation where there is 

very little concentration fluctuation and very large positive and 
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negative velocity fluctuations in the flow can be seen in Figure 
6. This data is taken in the recirculation zone where flow 
reversal occurs, but no large-scale motion is apparent. Data 
taken within the shear layer between the two jets is shown in 
Figure 7 , Here there is much scatter in the data and a strongly 
negative correlation. 

Once the detector was selected, a conditional sampling 
technique was used to determine the large-scale influence on the 
mass transport. Conditional sampling is a technique used to 
extract that portion of the data associated with the large-scale 
structures (as identified by the detector). One of the earliest 
uses of this experimental technique wass that of Kovasznay, et al 

(27) in 1970 on a turbulent boundary layer. Mathieu & Charnay 

(28) have reviewed this sampling method and much of the work done 
utilizing it. They point out that the difficulties of conditional 
sampling "are not due to the detection of signals", but are 
"consequences of the delicate choice of an adequate criterion". 
Blackwelder & Kaplan (24) used conditional averaging to examine 
the effect of coherent structures on the Reynolds stress for a 
turbulent boundary layer. The idea of conditional sampling or 
averaging, again from Mathieu & Charnay, is simply "the 
observation of a property only when some criterion is satisfied". 
This basic idea has been adapted to an analytical technique, 
rather than an experimental one, for this study. 

A computer program was written to perform conditional sampling 
on the simultaneous concentration-velocity data using 
concentration fluctuation as a trigger. A flow chart of this 
program is shown in Figure 8 and a program listing is included in 


Appendix A. This program made it possible to quantitatively 
separate the effect of the large-scale motion from the overall 
turbulent motion on the mass transport. Figures 9, 10, and 11 
show sample results from this program for the non-swirling flow. 
The section under DATA OUTPUT .... includes various statistical 
parameters calculated using all of the data at that point (these 
results are those reported in Ref. 11, 12). The lower section 
under CONDITIONAL SAMPLING RESULTS contains statistical parameters 
based on conditional sampling of the data according to 
concentration fluctuation. The information shown includes the 
concentration fluctuation range and the number of data samples in 
the range. Also included are mean velocity, relative mean 
velocity, root mean square velocity fluctuation, mass transport 
coefficient, and transport ratio, all calculated using only the 
data in the range. The relative mean velocity is defined as the 
mean velocity for a given concentration fluctuation range minus 
the mean velocity for all the data samples at that point. The 
transport ratio is defined as the mass transport coefficient: 

— i — i — 

f u 

*Vms u rms 


for the data in a given concentration fluctuation range divided by 
the overall mass transport coefficient. The local fraction of 
inner jet fluid is used for concentration (f) and f* = 1.0 by 
definition at the sudden expansion (Fig. 1). 

An example of data taken at a point on the inner edge of the 
shear layer is shown in Figure 9. The first interesting item is 
the negative overall mass transport coefficient. This combined 
with the fact that < 0 leads to the following: f a 
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3f 

3x 


>0 


which is inconsistent with normal gradient transport models 
already discussed. Moving to the conditional sampling results, it 
is seen that the majority of the large concentration fluctuations 
are negative. This is to be expected here since the inner jet 
fluid contains the dye (Fig. 1 ) and thus has f-1.0 . Any fluid 

swept into the region due to the large-scale motion will be outer 
jet fluid, with ^0. , and will result in a negative concentration 
fluctuation. An examination of the conditional sampling results 
shows two distinct flow characteristics among the data. The first 
is for the the data with I f | < • 2 , where the gradient transport 
model seems to fit. This data has a positive (or slightly 
negative) mass transport coefficient and very small relative mean 
velocity, indicating this well-mixed fluid is moving almost 
uniformly in the axial direction. The second flow characteristic 
is for the data with large concentration fluctuations (in this 
case |^-*2 ). This data is seen to have orders of magnitude, 

larger relative mean velocities and transport coefficients than 
the other data. the larger relative mean velocity indicates this 
fluid is not well-mixed with the rest of the flow and is moving at 
a considerably different, in this case higher, velocity than the 
mean flow due to the large-scale motion. The much larger mass 
transport coefficient indicates that even though the number of 
data samples is small, this data is responsible for the majority 
of the axial mass transport at this point. This data has a 
negative transport coefficient which, as already shown, is 
contrary to gradient transport. These data appear to be part of 
large-scale structures not well mixed with the overall flow and 
exhibit characteristics inconsistent with gradient transport. 
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Figure 10 shows data from a point well inside the shear layer 
under the strong influence of the large-scale motion. Once again 
the overall transport coefficient is negative, indicating that 
this is not a gradient transport situation. Since this location 
is inside the shear layer it is seen that there are large numbers 
of both positive and negative concentration fluctuations. There 
appear to be two distinct flow situations here also, but the 
large-scale motion seems to influence all of the data. Given that 
u'f = (Au) (Af) , notice in the conditional sampling results that 
when the concentration fluctuation (Af) is negative and the relative 
mean velocity (Au) is positive, u'f <0 an d hence a negative mass 
transport coefficient. The same result is observed when the 
concentration fluctuation is positive. The data within the 
large-scale structure can be separated by observing the large jump 
in relative mean velocity and transport coefficient. One again, a 
small number of data samples appears to be responsible for the 
jomority of the axial mass transport. If the data with f' <-0.2 
or f‘>.l is considered not part of the large scale, a weighted 
average of the large-scale data gives 36% of the samples 
accounting for 77% of the axial mass transport. The weighted 
average is calculated as the mean mass transport ratio for the 
large-scale data multiplied by the percent of the total number of 
samples that are included in the large-scale data. 

Finally, Figure 11 is an example of data containing no 
large-scale structure. There are no large concentration 
fluctuations and the relative mean velocitites are very small. 

This location is inside the recirculation zone in the outer region 
of the flow. 
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Precisely which concentration fluctuation ranges are included 
in the large-scale motion is an admittedly somewhat subjective 
decision. The conditional sampling program was run on axial and 
radial velocity concentration data at axial locations of 13, 51, 
102, 152, 203, 254, and 305mm for the non-swirling case. It was 
run at 13, 25, 51, 102, and 152mm for the swirling flow case since 
RobacX & Johnson (12) found that complete mixing occurred much 
more rapidly in this flow. In all cases, the data within the 
large-scale structure were separated out based on concentration 
fluctuation, as well as the clear jumps in relative mean velocity 
and mass transport coefficient. This technique, when applied to 
the Johnson & Bennett data set, resulted in a basically consistent 
set of results; it is thought that the scatter found in some of 
the results is due to the limited large-scale data available for 
the present investigation. Despite this, the hypothesis that 
large positive or negative concentration fluctuations can be used 
to identify large-scale motion appears to be valid. As stated 
earlier, the true test of a detector is the completeness of the 
pattern it yields. Once the concentration-velocity data for the 
large-scale structures were extracted, much information about the 
influence of these structures was available. The results, to be 
discussed, support this choice of detector and the resulting 
pattern is in fact well defined. 
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RESULTS 


As mentioned previously, the data sets used for the present 
investigation were acquired by Johnson & Bennett (11) for 
non-swirling flow and Roback & Johnson (12) for swirling flow. 

The results reported here are those associated with the 
large-scale structures. As discussed in the Analysis section, 
these results were obtained using conditional sampling techniques; 
the large-scale samples were detected using the concentration 
fluctuations as the criterion. 

Several parameters were examined using the output of the 
conditional sampling program. The complete set of results for the 
non-swirling flow at the axial station at 152mm is contained in 
Appendix B and consists of numerous sets of results for various 
r/R Q : axial velocity or radial velocity, each taken 

simultaneously with concentration. 

Using the procedures already described, several flow 

parameters were investigated as functions of r/R q for each axial 

location. The first parameter studied was the large-scale 

fraction defined as the ratio of samples with large concentration 

excursions compared to the total sample size. For the 

non-swirling flow, plots are included for the axial velocity 

(Figs. 12-19) and radial velocity (Figs. 20-25) data. In each 

case, the percentage appears to be maximized in the shear layer 

where the large scales are known to exist from flow visualization 

done previously (11). As already stated, the data acquisition 

rate was oriented toward repeatable overall statistics. Therefore 

the distribution between large-scale and non-large-scale data is 

believed to be representative of the actual flow. 
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It is seen that, with increasing downstream distance, the 
percentage of data in the large scale becomes a weaker function of 
radius, with the peak in the shear layer becoming less pronounced 
in the axial data. At 102mm (Fig. 14) there is a large peak in 
the shear layer at r/R Q = .2 where over 40% of the data is in the 
large scale. This drops off rapidly and by r/R q ? .4 is close to 
zero. At 152mm (Fig. 15) the curve is not as steep. From r/R Q = 

0 to r/ R q =.25, 30-40% of the data is in the large scale. The 
percentage approached 0 at r/R q - .5; but overall, more of the data 
is in the large-scale structure at this axial location. By 203mm 
(Fig. 16) the curve is still flatter. From r/R q = 0 to r/R Q ~ .2, 
50-60% of the data is in the large scale, not approaching zero 
until r/ R q =.6. Overall, a still larger percentage of the flow is 
part of the large-scale structure than in previous cases. 

Combining the axial velocity plots, a contour plot of percentage 
of data that is part of the large-scale structure was done (Fig. 
19). It is seen that the largest percentage is located within the 
shear layer and spreads radially, moving shownsteam as the shear 
layer spreads also. 

The radial velocity plots exhibit similar characteristics to 
the axial plots, with the peak spreading as the axial location 
increases. Figures 22 and 24 have curves that aren't quite so 
smooth as the others but the overall trend is similar. The 
percentage of data in the large scale must be the same for radial 
data as for axial; this is the case here for the most part. Any 
discrepancies can be attributed to the limited amount of 
large-scale data available. 
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Next, a calculation was made of the percentage of mass 
transport attributable to the large-scale portion of the flow. 

This was done for each radial location at a given axial station. 
This data, plotted as a function of t/Rq for each axial station, 
is shown in Figures 26 through 38 for the non-swirling flow. The 
axial mass transport percentage plots (Figs. 26-32) are typified 
by Figure 27, at z = 51mm, which shows that the percentage 
attributable to the large scale is a maximum within the shear 
layer between the jets. It is well in excess of 50% throughout 
the shear layer. 

Looking at the radial mass transport percentage plots at 13mm 
and 51mm (Figs. 33 and 34), they are similar to their axial 
counterparts (Figs. 26 and 27) although the percentage is higher 
for the radial data. Further downstream, these plots (Figs. 

35-38) were considerably more involved, again with higher peak 
percentages than for the corresponding axial ones (Figs. 28-31). 

Combining all of these results, it is possible to plot 
contours of percentage contribution to the mass transport due to 
the large-scale data. The axial mass transport plot (Fig. 39) 
shows that throughout the shear layer, 50-90% of the transport can 
be attributed to the large-scale motion. The increased complexity 
for radial mass transport became very evident when attempting to 
plot a similar contour plot for it. The limited locations with 
available data make it impossible to complete this figure. The 
existing data showed promise of a plot with similar features to 
the axial one, with somewhat higher percentages, but more data is 
needed to accurately complete this plot. 


21 


Looking at Figures 19 and 39, it is seen that there is 
qualitative agreement between the percentage of data that is part 
of the large-scale motion and the percentage contribution of the 
large scales to the mass transport. In the region where at least 
50% of the mass transport is attributed to the large-scale motion, 
at least 10% of the data is part of these large scales. In the 
peak region of 80-90% contribution, at least 20% of the data is 
part of the large-scale motion. These plots certainly 
substantiate the importance of large-scale motion in the mass 
transport of the flow. 

Figure 40 is a summary of the transport zones for axial mass 
transport. It is adapted from Figure 2 and includes the boundary 
of the region which Johnson & Bennett (10) found did not follow 
the traditional gradient transport model, as well as the boundary 
of the region for which the large-scale structure exerts a strong 
influence on the axial mass transport. The boundary of the 
large-scale region is based on a 50% contribution to the mass 
transport, although other values may be used by overlaying Figure 
2 with Figure 40. The boundaries of these two regions are 
observed to follow each other very closely. This leads to the 
conclusion that consideration of the large-scale structures is 
essential in accounting for the axial mass transport in this 
region where the gradient transport model is inadequate. Of 
course, the presence of the large-scale structures means the flow 
includes multiple scales; gradient transport should not be 
expected to apply in such regions. 
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For the swirling flow, plots of large-scale fraction (Figs. 
41-43) and mass transport percentage (Figs. 44-46) are included 
for axial velocity at axial stations of 13, 25, and 51mm. As 
Roback & Johnson (12) found, mixing for the swirling flow case 
occurs in approximately one third the downstream distance that it 
takes for the non-swirling flow. For this reason, large-scale 
structures were found to be influential in a much smaller region 
than for the non-swirling flow. 

At 13mm (Fig. 41) the large-scale fraction peaks at 
approximately 30% of the data samples while at 25mm (Fig. 42) the 
peak is closer to 40%. In both cases the maximum occurs within 
the shear layer between the inner and annular jets as expected. 

By 51mm (Fig 43) the large-scale fraction is peaking at over 60% 
but the peak has shifted to the centerline. This is apparently 
due to the unsteadiness associated with the leading edge of the 
recirculation zone which occurs for the swirling flow and is 
consistent with the findings of Roback & Johnson (12). As 
expected, these percentages are somewhat higher at the same axial 
locations for the swirling flow than for the non-swirling flow. 

The largest difference between swirling and non-swirling flow 
occurs at the 51mm location. The high values there are almost 
certainly associated with the unsteady motion of the recirculation 
zone noted earlier. Flow visualization (12) at this point clearly 
shows the very large structures present near the centerline. 

Also, as seen in the non-swirling flow, the overall percentage of 
data within the large-scale motion increases with downstream 
distance . 
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Swirling flow data were analyzed for large-scale structures at 
the 102 and 152mm locations. Though some indication of large 
structures was found, it represented, at most, 3% of the total 
samples at any one location. These results are not included as 
the statistics contain significant uncertainty due to the small 
sample size. We would agree with Roback & Johnson that 
large-scale structures exist at 102nun at least; but if the sample 
set is representative of the large-scale fraction, the large-scale 
structures are of no significance at this point. 

The swirling flow mass transport percentage plots (Figs. 

44-46) also exhibit some similarities to the non-swirling case. 

In each case the percentage is a maximum in the shear layer 
between the jets and the region of large-scale influence spreads 
radially with downstream distance. As already stated, the 
swirling flow mixes very rapidly and these plots reflect this. 

The percentages are high, like the non-swirl flow, but spread into 
the centerline much sooner axially than for the non-swirl flow. 

A result of the rapid mixing for the swirling flow is that 
changes from station to station are much larger than for the 
non-swirl case. This makes it difficult to draw contours. As 
with the non-swirling flow, the influence of the large-scale 
motion is large and also consistent throughout the flow field. 

An effort was undertaken to utilize the two-component velocity 
data to look at the large-scale contribution to the momentum 
transport. It was necessary to select a new detector to separate 
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the large-scale data since there is no concentration data 
available in this case. A modified conditional sampling program 
was written to try to find such a trigger. A flow chart is 
included in Figure 47 and a program listing is in Appendix A. 

This program was used on the concentration-velocity data 
previously sampled for the non-swirling flow using concentration 
fluctuation as a trigger. Calculations were made of a conditional 
mean velocity and rms velocity fluctuation based on the same 
concentration range already selected. These were then compared to 
the overall mean and rms for the data. A typical output from this 
program is shown (Fig. 48) together with a plot of the data (Fig. 
49). This figure shows the velocity distribution, in this case, a 
radial velocity data set. The overall data exhibits the nearly 
normal distribution. The large-scale data is more skewed, with a 
significantly different mean velocity and somewhat different rms 
fluctuation velocity than the overall. It had been hoped that 
radial velocity or velocity fluctuation might be used as a 
detector for the large-scale structure in the two-component 
velocity data. However, based on the results typified by this 
figure, there is no way to use either to separate the large-scale 
data. Use of axial velocity as a trigger was ruled out also when 
plots of this type showed similar results. 

Plots of mean velocity and rms velocity fluctuation as 
functions of r/Rg were done for simultaneous non-swirling 
concentration-velocity (axial and radial) data at various axial 
stations. The overall and large-scale values were plotted 
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together for comparison purposes. Figure 50 is a typical plot for 
axial velocity. The mean velocity for the overall data and that 
for the large-scale data have very small differences, indicating 
that the coherent structures are moving at approximately the same 
axial velocity as the rest of the flow. This result for free 
shear layers is consistent with the results of Kovasznay et al 
(27) and several others (reviewed by Cantwell 18) who found large 
structure convection velocities in wall-bounded flows between 80 
and 90 percent of the free stream mean velocity. Roshko (17) 
noted that for plane-mixing layers the large-scale vortices 
convect at approximately the average velocity of the two streams, 
with the velocity and density ratios having some influence. 

Plots of mean radial velocity vs. r/Rp are included for axial 
stations of 51, 102, 152, and 203mm. Figure 51 shows a negligible 
difference between the large-scale and overall data at 51mm. At 
102mm, the differences are very significant (Fig. 52). For 
0_<r/RQ<.2, the large-scale structures have a larger negative 
radial velocity than the mean flow. From . 2_< r /R q£. 4 , the 
large-scale structures have a larger positive radial velocity than 
the mean flow, with the difference increasing with r/R Q . In this 
region, the shear layer is spreading radially, assuming that the 
large-scale structures identify the shear region. 

Further downstream, at 152mm (Fig. 53), there is a negligible 
difference between large-scale motion and overall data for 
t/Rq<. 2. This difference is significant, however, for 
. 275^r/RQ^. 475 and is approximately constant in this range. Again 
the coherent structures have a significantly larger radial 
velocity than indicated by the mean velocity of all the data. 
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At 203mm (Fig. 54), there are very small differences for 
0 <t/Rq<^ 3. Large differences are found between the large-scale 
flow and overall for . 3<r/RQ^.6. The radial mean velocity appears 
to be a linearly increasing function of r/Rg for the large-scale 
and overall data, although the large-scale increase at a faster 
rate. 

If the extent of the shear layer is assumed to correspond to 
the region where large-scale structures are found, the spreading 
rate of the shear layer for the non-swirling flow can be examined 
using Figures 12 through 17. Table 1 shows the radial extent of 
the large-scale structures for each axial station taken directly 
from the figures. Data were not taken at the same radial 
locations for all axial locations; therefore, comparisons of the 
exact same radial locations are not always possible. The data 
base is complete enough to observe that the spreading rate of the 
shear layer, or region where large-scale structures exist, 
increases approximately linearly with increasing axial location. 
This is consistent with other work on free shear layers (29) and 
consistent with the positive radial velocity found by Johnson & 
Bennett (11). These results emphasize the need for multiple scale 
modeling (15) since the large scales and overall flow have 
different radial convection velocities. 

Figures 55, 56, and 57 are plots of the large-scale relative 

mean velocity as a function of concentration fluctuation, at axial 

locations of 51, 102, and 152mm respectively, for the non-swirling 

flow. The large-scale relative mean velocity is defined, at each 

axial and radial location, as the mean velocity for a given 

concentration fluctuation range minus the mean velocity for all of 

the large-scale data at the location. These plots provide 
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additional support for the selection of the large-scale structures 
in that the large-scale relative mean velocity is approximately 
zero, in each figure, for concentration fluctuation of zero. 

At any location there is some fluctuation in the overall 
curve, attributed to eddy variations. Each curve has limited 
scatter though and the vortex like structures are quite well 
defined in that a given relative mean velocity suggests a 
preferred concentration fluctuation value. Since the structures 
are well defined, the lack of correlation reported by Johnson & 
Bennett (11) must be due to the random arrival times of the 
structures. This is also suggested by the 3-D calculations of 
Riley & Metcalfe (30) and leads to the conclusion that the present 
form of the structures is "3-D like" for the large-scales. 
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CONCLUSION 


The initial or upstream development of shear flows is crucial 
to the performance of many aerodynamic components (i.e., 
combustors). Comparisons of experimental results with computer 
predictions show that present models are not sufficient for such 
flows. Historically, flow visualization has shown that 
large-scale structures are present within this upstream, 
non-equilibrium region. It has been generally hypothesized that 
any lack of agreement between prediction and experiment can be 
attributed to these structures. In general, no adequate model 
(including them) has been included in the numerical computer 
codes, primarily since little was known about them. 

The present investigation set out to analyze the physical 
characteristics of large-scale structures in free shear layers and 
to evaluate their importance to turbulent transport (including the 
first ever results for mass transport) . Utilizing previously 
acquired data, the conclusions are limited in that the data set 
was not optimized for this use. Regardless, several significant 
conclusions are possible; they include: 

(1) The large concentration fluctuations are indicative of 
the large structures. First, the sign of the transport is 
predictable with a model consistent with these structures. 
Secondly, the region with significant large concentration 
fluctuations is consistent with the large-scale region noted from 
flow visualization. 
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(2) The large structures, where present, account for most of 
the mass transport; the obvious implication of multiple scales is 
consistent with the lack of agreement with a gradient transport 
model . 

(3) The large structures are found to be convected axially 
with the overall mean velocity but radially at a faster speed; 
this further strengthens the need for a multiple scale model. 

(4) Any investigation of the influence of large-scale 
structures on momentum transport in three-dimensional free shear 
layers will most probably require concentration detection. 

(5) The effects of large-scale motion for the swirling flow 
case or limited, at least for the data available, to the region 
just upstream of the recirculation zone and to the shear region 
between the jets. 

The present investigation therefore has documented the 
importance of the large-scale structures and the need for them to 
be modelled if computer predictions are to be accurate. It is 
strongly suggested that the development of such a model is of the 
highest priority. Experimentally, this suggests several projects 
for further investigation. They include: 

(1) A detailed analysis of the flow development with the 
large structures. What is the structure in physical space? Whay 
type of mixing occurs within them (is it a two-step process as 
suggested for two-dimensional shear layers)? 

(2) A detailed analysis of momentum transport associated with 
the large structures. Does if significantly affect the results 
as would be predicted from the present mass transport results and 
from the results for other shear flows? 
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(3) An investigation of the influence of initial conditions 
on large-scale development. Certainly, such initial conditions 
are required for any computer simulation. What effects do changes 
in initial conditions have and what are the implications for 
accurate models? 

(4) An investigation of two-dimensional vs. 

Three-dimensional shear layers. Both the published experimental 
results and the recent computer simulations strongly suggest that 
such a delineation between two-dimensional and three-dimensional 
flows is not straightforward. What are the implications of such 
assumptions on the development of an accurate model? 
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TABLE 1 
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Fig. 1 


SHEAR REGIONS OF COAXIAL JETS 




reprinted from NASA CR-168252 


33 


Fig. 2 
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END OF RUN 61 -POINT 
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END OF RUN 61 -POINT 

terminated: stop 
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END OF RUN 52 -POINT 14 

terminated: stop 
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Fig. 14 
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OF POOR QUALITY 

Fig. 35 


RADIAL COMPONENT AT Ze102 MM 
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MASS TRANSPORT PERCENTAGE 


Fig. 37 
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RADIAL COMPONENT AT Z-203 MM 
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AXIAL MASS TRANSPORT 
LARGE SCALE CONTRIBUTION 



AXIAL POSITION ,2 (mm) 


TRANSPORT ZONES SUMMARY 


Fig. 40 
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AXIAL POSITION , Z (mm) 


PERCENTAGE OF DATA IN LARGE SCALE 


Fig. 41 
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PERCENTAGE OF DATA IN LARGE SCALE 


Fig. 42 



RADIUS RATIO , (r/Ro) 
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MASS TRANSPORT PERCENTAGE 


Fig. 45 
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Fig. 46 
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AXIAL COMPONENT AT Z=5 1 MM 
SWIRLING FLOW 
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Fig. 47 



/ READ l UFO PROM- HEADER ( x , r / R o , • t c . ) 
READ ALL DATA (velocity A concentration! 



INPUT CONCENTRATION FLUCTUATION 
RANOe FOR LARQE-SCALE OATA 


SAMPLE ALL OATA: CALCULATE NUMBER OF LAROE'SCALE 
SAMPLES, MEAN A RMS VELOCITY FON LA ROE-SCALE 

OATA: LOCATE MIN A MAX VELOCITY 


DIVIDE OATA INTO EQUAL INTERVALS BASSO ON VELOCITY 
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Fig. 


*ROf 31 .15 

SUCCESSFUL TEnP » RESTORE RUH31P13 *S (TPl5) 
CANCELLED: 30 na«E RrtOS unknown 


data OUTPUT FOR RUN SI FOINT 13 


RADIAL VELOCITY vs CONCENTRATION 
*•132 Art ANO ft/ RO- 0.324 

SELECT HlNlrtUN POSITIVE AND NEGATIVE CONCENTRATION FLUCTUATION 
TO IE INCLUDED IN LARGE SCALE STRUCTURE USING F4.1 FORHATl 
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0.2 

CFNEG- 

• 1.0 


nuhDER of SAMPLES i ORIGINAL* 


999 GOOD DATA-984 
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RADIAL VELOCITY DISTRIBUTION 
Z= 152 mm r/Ro-.324 
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RADIUS RATIO (r/Ro) 
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VELOCITY (m/s) 


Fig. 51 


RADIAL MEAN VELOCITY Z * 6 1 mm 


LARGE SCALE 

OVERALL 
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VELOCITY (m/s) 


Fig. 52 


RADIAL MEAN VELOCITY Z=102 mm 


LARGE SCALE 
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VELOCITY (m/3) 


Fig. 53 


RADIAL MEAN VELOCITY Z=152 mm 


LARGE SCALE 

OVERALL 

r/BoiO ° 
r/Ro< 0 ■ 

r/Ro>0 O 
r / R o<0 • 
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VELOCITY (m/s) 


Fig. 54 


RADIAL MEAN VELOCITY Z=203 mm 



RADIUS RATIO (r/Ro) 
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NO N - SWIRLING FLOW 
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HOLD * CV 


APPENDIX A 


OOOOIOO 

0000200 

0000300 

0000400 

0000500 

0000600 

0000700 

0000733 

0000770 

0000790 

0000800 

0000900 

0001000 

0001100 

0001130 

0001200 

0001300 

0001400 

0001500 

0001600 

0001700 

0001800 

0001900 

0002000 

0002100 

0002130 

0002200 

0002300 

0002400 

0002500 

0002600 

0002700 

0002800 

0002900 

0003000 

0003035 

0003070 

0003100 

0003200 

0003300 

0003400 

0003500 

0003600 

0003700 

0003800 

0003900 

0004000 

0004100 

0004200 

0004300 

0004400 

0004500 

0004600 

0004700 

0004800 

0005000 

0005100 

0005200 

0003300 

0005400 

0005500 

0005600 

0005650 

0005700 

0005800 

0005900 

0006000 

0006100 

0006200 

0006300 

0006400 

0006500 

0006600 

0006700 

0006800 

0006900 

0007000 

0007100 

0007200 

0007300 


DI HEMS I ON V< 1000) »C( 1000) »VN< 1000) 

REAL*8 VEL 
REAL NEVBAR*NEURMS 
INTEGER RUN* POINT tZ 
VBAR-O. 

VRMS-O. 

N4*0 

THIS PROGRAM PERFORMS CONDITIONAL SAMPLING ON RUNS OF CONCENTRATION- 
VELOCITY DATA( METRIC UNITS)* 

READ INFO FROM HEADER 
READ< 11 . 10) RUN * POINT 
10 FORMAT<26XtI2»10X»I2) 

READ ( 1 1 * 1 5 >N0 
13 FORMAT < 4 1 X * I 3 > 

NFIRST-NO 
READ< 1 1 »20)Z.RR0 
20 FORMAT ( 27X # I 3 » 1 2X » F6 . 3 ) 

READ ( 1 1 >25) 

READ( 11 >25) 

25 FORMAT <5X ) 

READ(1 1 » 30) VEL 
30 FORMAT ( 34X » 1A8) 

DO 40 N*1 * 4 
READ< 1 1 >25) 

40 CONTINUE 

READ ALL DATA) CALCULATE INITIAL MEAN 1 RMS 
DO SO I * 1 * NO 
REA0( 1 1 * 45) V< I ) »CU) 

45 FORMAT (33XfF6*3»7X»F6*3) 

VBAR-VBARFV< I ) 

VRMS*VRMS4V( I ) **2 
50 CONTINUE 

VBAR«VBAR/FLOAT(NO) 

VRMS*VRMS/FLOAT < N0>-VBAR**2 
VRMS«SQRT(VRMS) 

THROW OUT DATA FARTHER THAN 3*VRMS AWAY FROM VRAR t CALCULATE 
NEW MEAN t RMS WITH DATA LEFT OVCR . 

DO 200 N* 1 > 5 

IF( <FL0AT<N4)/FL0AT<N0) ) .GT *0.99)00 TO 700 
IF < N . GT , 1 )N0*=N4 
NEWB AR -*0 • 

NEURMS-O* 

N4*0 

DO 100 1*1 # NO 

IF <ABS(V< I )-VBAR) *GT . <3 . 0 * VRMS ) )G0 TO 100 
N4 ® N 4 + 1 

NEWBAR*NEWBARiV( I ) 

NEWRMS=N£URMS+V< I )**2 
V(N4)-V< I ) 

C ( N4 ) *C ( I ) 

100 CONTINUE 

VBAR*NEWBAR/FL0AT(N4) 

VRMS*NEURMS/FLOAT <N4)-VBAR**2 
VRMS*SQRT ( VRMS) 

200 CONTINUE 
Vl*0. 

V2*0 . 

U3 = 0 • 

U4 *0 . 

01 * 0 . 

02 * 0 . 

G1 *0. 

CALCULATE ASSORTED STATISTICAL PARAMETERS < MEAN > RMS » MOMENTS* ETC . > 

DO 500 I = 1 * N4 
V1»V1+V< I ) 

V2*V2tV< I ) **2 
U3-U3+V ( 1)4*3 
U4 =U4 ♦ V ( I ) * *4 
01=Q1+C< I ) 

02*Q24C < I > **2 
G 1 *G 1 ♦ V ( I ) *C ( I > 

500 CONTINUE 

V 1 *V 1 /FLOAT ( N4 ) 

U4*U4/FL0AT(N4 ) -4 . *V1 *U3/FL0AT<N4>+6. *V1**2*V2/FL0AT <N4>-3. tVl **4 
U3*U3/FLQAT<N4 >-3 * *V1 *V2/FL0AT (N4 >42 . *V1**3 
V2-V2/FL0AT <N4 > -VJ **2 
U2-S0RT ( V2) 

01 *01 /FLOAT l N4 ) 

O2-Q2/FL0AT (N4 ) -01**2 
02*SPRT < 02 > 
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0007400 
0007300 
0007000 
0007700 
0007800 
0007900 
0008000 
0008100 
0008200 
0008300 
0008400 
0008300 
0008330 
0008300 
0008390 
0008393 
0008000 
0008700 
0008800 
0008900 
0009000 
0009100 
0009200 
0009300 
0009400 
0009300 
0009600 
0009700 
0009800 
0009900 
0010000 
0010100 
0010200 
0010300 
0010400 
0010300 
0010600 
0010700 
0010800 
0010900 
0011000 
0011100 
0011130 
0011200 
0011300 
0011400 
0011500 
0011600 
0011700 
0011800 
0011900 
0012000 
0012100 
0012200 
0012300 
0012400 
0012500 
0012600 
0012700 
0012800 
0012900 
0013000 
0013100 
0013200 
0013300 
0013400 
0013300 
0013600 
0013700 
0013800 
0014400 
0014S00 C 
0014600 
0014700 C 
0014800 
0014900 
0015000 C 
0013100 
0013200 C 
0015300 


01 *01 /FLOAT CN4 ) -Q1 tVl 

N6-0 

P8-0. 

R9-0. 

N7-0 

R8-0. 

809 URITE<7*810) 

810 FORMATC/) 

UR I TEC 7 * 8 1 5 )RUN »P0INT 

813 F0RMATC3X* 'DATA OUTPUT FOR RUN'*X3*' P0INT**13) 

WRIT£<7*820) 

820 FORMAT < /) 

URITEC 7*622) VEL 

822 FORMAT (1X*1A7*' VELOCITY VS CONCENTRATION* > 

URITEC 7*823 )Z*RRO 

823 FORMAT ClX*'Z**tI3#' MM AND R/R0***F6.3) 

URITEC 7*823) NFIRST *N4 v 

823 FORMAT (IX* 'NO*' 14* * AND N4*'*13) 

UR I TEC 7 v 830) VI 

830 FORMAT C 1X» * VBAR* ' * F9»4* * MPS*) 

URITEC7f 840>U2 

840 FORMAT C lXf ' VRMS* '*F9*4»' MPS') 

MRITE(7»850>U3 4 

830 FORMAT (IXt'THIRD MOMENT 0* TURBULENCE** *E12*S*' MPS443') 

R3-U3/U2B83 

WRITE(7»860)R3 

860 FORMAT < IX * 'THIRD CORRELATION COEFFICIENT** *F9.4> 

URXTEC7* 8701U4 

870 FORMAT C 1X» 'FOURTH MOMENT OF TURBUl ENCF-' *E12.5* * MPSBB4 * ) 

R4*U4/U2tS4 
URITEC 7*880 )R4 

880 FORMATC IX* 'FOURTH CORRELATION COEFFICIENT** *F9.4) 

URITEC 7*881 >Q1 

881 FORMAT (IX* *C8AR* ' *F8 . 3 * ' X ' ) 

URITEC 7*882)02 

882 FORMAT C IX* *CRMS* '*F8.3*'X') 

UR I TEC 7*683)01 

883 FORMAT C IX* 'CPVPBAR* **F10.6) 

TT-G1/Q2/U2 

URXTEC 7* 884 ) TT 

884 FORMAT C IX* 'OVERALL TRANSPORT COEFFICIENT*' #F10.6> 

URITEC7*810) 

CONDITIONAL SAMPLING SECTION BEGINS HERE I 
URITEC7*890) 

890 F0RMATC29X* 'CONDITIONAL SAMPLING RESULTS*) 

URXTEC7*900) 

900 FORMAT C / ) 

URXT£C7*910> 

910 FORMAT (IX* 'CONCENTRATION' *4X. 'NUMBER OF '. 13X* 'RELATIVE '* 19X* ' TRANSPORT ' *6X *' TRANSP0R1 ' ) 
URITEC7*920> 

920 FORMAT C2X* 'FLUCTUATION ' * 5X* 'OCCURANCES ' * 4X* 'MEAN' *8X* 'MEAN* * 

! 10X* 'RMS' 7X» 'COEFFICIENT' »7X* 'RATIO') 

DO 973 K*1 * 20 
Al — 1.0+<K-l)t0.1 
A2-A19.0999 

Ul-O. 4 

NCOND-O 

V3*0 • 

V4*0* 

VR*0* 

DO 923 N*1 *N4 

TEST TO MAKE SURE C' IS IN THE DESIRED RANOE • y 

IF(V<N).EQ.-100. > GO TO 925 
IFCCCCN)~Q1)*LT*A1) GO TO 925 
IF(<CCN)-Q1).GT.CA2) ) GO TO 925 
V3-V3+VCN) 

V4-V4+VCN)M2 

NCOND-NCONDil 

Ul«Ut4CC(N)-0t)t(V(N)-Vl) 

923 CONTINUE 

927 IFCNCOND.EQ.O) GO TO 950 
03 IS CONDITIONAL MEAN VELOCITY. 

V3*V3/FL0AT < NCOND > 

V4 IS CONDITIONAL RMS VELOCITY FLUCTUATION. 

V4*V4 /FLOAT C NCOND ) -VJM2 
V4-S0RTCV4) 

U1 IS CONDITIONAL MASS TRANSPORT COEFFICIENT. 

U 1-U1 /FLOAT (NCOND) /U2/Q2 

VR IS RELATIVE MEAN VELOC I TT ( COND I T I ONAL MEAN-OVERALL MEAN). 

VR-V3-VI 
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0015400 

0013500 

001S600 

0013700 

0013B00 

0015900 

0016000 

0016100 

0016200 

0016300 

0016400 

0016300 

0016600 

0016700 

0016800 

0016900 

0017000 


Y IS INC HASS TRANSPORT COEFFICIENT RAT 10! CONDI T TONAL /OVERALL) • 

Y-Wl/TT 

WRITE! 7»930)A1 »A2»NC0ND»V3*VR»V4 » W1 r Y 

930 FORMAT! IX »F4 • 1 » ' - ' • ► /,4*AX»IJ»3X»F8.4»4X»F9.3f4X»t'9.3»4X»fl0.5»5X»Fl0.3> 
00 TO 974 

950 WRITE < 7 » 960 > A 1 * A 2 • NCOND 

960 FORMAT !lXtF4.1»' - ' * F 7 . 4 » 6X 1 1 3 ) 

X IS THE CENTER OF THE C' INTERVAL • 

974 X-A19.05 

92B IS THE FORMAT STATEMENT US1D TOR WRITING PLOT FILES. 

928 FORMAT !F5.2»F10«3) 

973 CONTINUE 

WRITE! 7»H20> 

963 WRITE! 7 » 970) RUN# POINT 

970 FORMAT ! IX» 'END OF RUN'.I3»' -POINT' *13) 

STOP 

END 
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OOOOIOO DIMENSION Nl C 1000 )» N2 < 1000 ) 

0000200 DIMENSION V!1 000 > t C <1000 U VN ( 1 000 ) 

0000300 REAL t8 VEL 

0000400 REAL NEW BAR • HE UR MS 

OOOOSOO INTEGER RUNtPOINTrZ 

0000600 VBAR-O. 

0000700 VRMS-O, 

0000800 N4»0 

0000900 C THIS PROGRAM PERFORMS CONDITIONAL SAMPLING ON RUNT* OF CONCENTRATION- 
0001000 C VELOCITY OA TA ! METRIC UNITS) AND CALCULATES PROBABILITY DENSITY DISTRIBUTIONS 
0001100 C READ INFO FROM HEADER 


0001200 


READ! 1 1 » 10) RUN* POINT 


0001300 

10 

FORMAT (26Xr!2»10XvI2> 


0001400 


READ! 1 1 » 1S)N0 


0001500 

15 

FORMAT !41X*I3) 

5> 

0001600 


NFIRST-NO 


0001700 


READ! lit 20)2 rRRO 


0001800 

20 

FORMAT !27X»I3tl2X»F6. 3) 


0001900 


READ! 11 *25) 


0002000 


READ! 11*25) 

t 

0002100 

25 

F0RMATC5X) 


0002200 


READ! lit 30) VEL 


0002300 

30 

FORMAT! 34Xf IA8) 


0002400 


DO 40 N-l»4 


0002500 


READ! 11 r 25) 


0002600 

40 

CONTINUE 



0002700 C READ ALL DATA! CALCULATE INITIAL MEAN I RMS 
0002800 DO 50 I«lrNO 

0002900 READ( 1 1 » 45 ) V ( I ) »C< I ) 

0003000 45 F0RMAT!33XrF6.3t 7X*)-4.3> 

0003100 VBAR-VBAR9V! 1 ) 

0003200 VRMS-VRMS+V! I >tt2 

0003300 SO CONTINUE 

0003400 VBAR«VBAR/FL0AT!NO) 

0003500 VRMS«VRMS/FL0AT!N0)-VBARtt2 

0003600 VRMS- SORT ( VRMS > 

0003700 C THROW OUT DATA FARTHER THAN 3BVRMS AWAY FROM VBAK S CALCULATE 
0003800 C NEW MEAN t RMS WITH DATA LEFT OVER. 

0003900 
0004000 
0004100 
0004200 
0004300 
0004400 
0004500 
0004600 
0004700 
0004800 
0004900 
0005000 
0005100 
0005200 
0005300 
0005400 
0005500 
0005600 
0005700 

oooseoo 

0005900 
0006200 
0006500 
0006600 
0006700 
0006800 
0006900 
0007000 
0007100 
0007200 
0007300 
0007400 
0007500 
0007600 
0007700 
0007800 
0007900 
0008000 
0008100 
0008200 
0008300 
0008400 
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DO 200 N-1'5 

IF! (FLOAT (N4 ) /FLOAT < NO) ) • G l •0.99)00 TO 200 

IF<N'GT.1)N0»N4 

NEUBAR-O. 

NEWRMS-O. 

N4*0 

DO 100 I * 1 * NO 

IF<ABS<V<I )-VBAR) «GT • ! 3*0tVRMS> )G0 TO 100 
N4*N44 1 

NEWBAR-NEWBAR4V<I> 

NEWRMS-NEWRMSf V < I > 1 1 2 
V(N4 > *V< I ) 

C<N4 ) »C< I ) 

100 CONTINUE 

VBAR«N£WBAR/FLOA T < N4 > 

VRMS*NEWRMS/ FLOAT <N4)-VBARtt2 
VRMS-SORT(VRMS) 

200 CONTINUE 

00 500 I» l * N4 
Vl«VltV!I> 

V2»V2+V!I)t«2 

Q1-Q14C<I) 

500 CONTINUE 

VI • VI /FLOAT < N4 > 

V2-V2/FL0AT < N4 > -VI 1 12 
U2-SQRT<V2) 

Q1-Q1/FL0AT<N4) 

809 WRITE(7»810> 

810 FORMAT < / ) 

WRITE! 7 t 61 5) RUN » POINT 

815 FORMAT <5X i 'DATA OUTPUT FUR RUN'rIJr' P0INT'#13> 

URITE<7»820) 

820 FORMAT!/) 

WRITE ( 7*822) VEL 

822 FORMAT! IX. 1A7.' VELOCITY VS CONCENTRATION') 

URITE<7t823)Z»RR0 

823 F0RMAT<1X» 'Z-' »I3. ' MM AND R/RO«* ' * F 6 . 3 ) 

WRI TE< 7» 25 > 

WRI TE < 7 » 980 > 

980 F0RMAT!3X» 'SELECT MINIMUM POSITIVE AND NEGATIVE CONCENTRATION FLUCTUATION') 
WRITE<7»901 ) 

981 FORMAT < 1 X» 'TO BE INCLUDED IN l ARGE SCALE STRUCTURE USING F4.1 FORMAT I ' ) 


0008300 

0008600 

0008700 

0008600 

0008900 

0009000 

0009100 

0009200 

0009300 

0009400 

0009300 

0009600 

0009700 

0009800 

0009900 

0010000 

0010100 

0010200 

0010300 

0010400 

0010300 

0010600 

0010700 

0010800 

0010900 

0011000 

0011100 

0011200 

0011300 

0011400 

0011500 

0011600 

0011700 

0011800 

0011900 

0012000 

0012100 

0012200 

0012300 

0012400 

0012300 

0012330 

0012600 

0012700 

0012800 

0012900 

0013000 

0013010 

0013020 

0013030 

0013040 

0013030 

0013060 

0013070 

0013100 

0013200 


WRITE(7*982> 

982 FORMAT (2X* 'CPPOS* ' ) 

READ! 7 *983) CPPOS 

983 F0RMAT(F4.1) 

UR1TEC7* 984 ) 

984 FORMATdX. 'CPNEG-') 

R£AD(7*983)CPMEG 

M5»0 

UBAR-O. 

URMS-O. 

DO 999 1*1 fN4 
Ml ( I 1*0 
M2< I )*0 
999 CONTINUE 

VMIN*VBAR-3 • OtVRMS 
VMAX*VBAR43 • OtVRMS 
DO 1000 N-1.N4 

IF((!C(N)-QI)#LT . CPPOS) * AND* C (C(N)-Ol) . GE* CPNEG) 1 CO TO 950 
N3-M341 

U8AR*U8AR4V(N) 

URMS*URMS4V(N)t«2 

I*XNT( < ( V(N)-VMIN) / ( VMAX-VNIN) )9S0)41 
iFd.EO. 51)1*50 

C N1 IS THE NUMBER OF CONDI1IOMAL SAMPLES 
Nl(I)*N 1(1)41 

950 1»INT < ( <V(N)-VMIN)/(VMAX-VM1N>)»50>41 
IF ( I *CQ»31 ) 1*30 

N2 IS THE NUMBER OF TOTAL SAMPLES 
N2( I ) *N2( 1)41 
1000 CONTINUE 

UBAR k URNS ARE CONDITIONAL (LARGE SCALE) MEAN SRMS VELOCITIES 
UBAR-UBAR/FL0AT!N5> 

URMS-(URHS/FL0AT(N3)-UBAR»«2)**.S 
WRITE! 7*25) 

WRITE ( 7* 1010)NF1RST »N4 r N5 

1010 FORMATdX* 'NUMBER OF SAMPLESI ORIGINAL* ' * 13 * ' GOOD DATA-'*I3.' LARGE SCALE-'. 13) 
WRITE (7*25) 

WRITE< 7» 1025) 

1023 F0RMAT!25X* 'NUMBER OF SAMPLES') 

WRITE! 7* 1030 > 

1050 FORMAT <3Xt ' I ' *5X * ' VMID(M/S) ' *3X* 'L.S. ' * 5X* ' OVERALL ' ) 

WRITE! 7*25) 

DO 1100 1*1*30 

VHID*VMIN44UHAX-VMIN)/50. *<!-.£> 

WRITE ( 7* 1073) I*VMXD*N1!I)»N2< X) 

1073 FORMAT (4X*12*6X*F6«3*7X*X3*7X»13> 

1100 CONTINUE 

WRITE<7*25> 

WRITE! 7* 1300 >VBAR* UBAR 

1300 FORMATdX* 'MEAN VELOCITY ! M/S ) I OVERALL* ' * F6 • 3 * ' I LARGE SCALE* ' *F6 . 3) 

WRITE (7» 1600) VRMSfURMS 

1600 FORMATdX* 'RMS VELOCITY!M/S) I OVERALL* ' * F 6. 3 • ' I LARGE SCALE-' *F6.3) 

WRITE! 7* 1700)01 -CPPOS tCPNEG 

1700 FORMAT (IX » ' CBAR*' *F6 • 3* ' CPPOS*' *F6 . 3 * ' CPHEO-' »F6 . 3) 

STOP 

END 
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SAMPLE 63.17 

SUCCESSFUL < TEMP ) RESTORE RUN63P17 AS <TP17> 

cancelled: ddname rmds unknown 
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END OF RUN 63 -POINT 17 

terminated: stop 
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SAMPLE 63.19 

SUCCESSFUL < TEMP ) RESTORE RUN63M9 AS (IF 19) 
CANCELLED? DONANE RMDS UNKNOWN 
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END OF RUN 63 -POINT 19 

terminated; stop 


SAMPLE 63.20 

SUCCESSFUL (TEMP) RESTORE RUN63P20 AS (TF20) 
CAHCetLLD: UUHAHZ fiHb S UNKNOWN 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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END OF RUN 63 -POINT 20 

terminated: stop 


SAMPLE 6irJl 

SUCCESSfUL (TEMP) RESTORE RUN63T21 

cancelled: ddname rmds unknown 


►“ 

cc 



MAN 

n o ci 

o 

► 

r* -o 

o 

O 

♦ 

dm n 

© O' CD 

w >o 

fs 

n rs 

Cc 


O 

r* r» o 

<r 

— CD 

4T 

© d 

in 

>- 

• 

♦ ♦ * 

♦ * • 

• ♦ 

• 

• • 

z 



0 

0 

t 

t 

0 

0 

c* n 

<r 

’■O o 

<t 

£ 

*- 

£ 








»- 

oc 

£ 

Ul 

d 

rs ® 

© 

© -* o* 

0* 

Kj 

O' 

rs *■> 

o 


rs 

-• +■ 

© 

rv n o 

o» 

'O 

o 

r. ro 

a. 

u 

'O 

cy n 

♦ 

^ © -4 

m 


K> 

rs -i 

cn 


rs 

n cd 


IT- ^ « 



SJ 

rs w* 

z 

u. 

© 

/S C 4 

© 

O d rs 


W 

rt 

rt 

<z 

u. 

» 




♦ 

• 



£ 

ui 

*4 

0 

0 

© 

0 

0 

0 

v4 

w 4 

d 

n n 


o 

a 

' 

i 


t . i 

* 

1 


i > 


^-vaD(^ocDon«H^oro 
GiriN^(L'Or>,HfflNno 
t n 0'nnn+-nnnc4f*.ncD 

r r<ON»C'ON'0«<^H 

OC f<MH«MN***<*<***<*< 


(A 


OOOOOOOOOOOO 


M 

4. 


<s> 

<t 









© 

z> 

flsrinno'ON’O’OrioA 




■* 


z 

~ z 

O' — o 




K» 

m 



►- <z 

00 < o 

O-OONSC^-flA 



w 

m 


-1 

<X Ul 

P 4 ^ 




f* 

c n 


Cl 

-1 £ 





tn 

£ 


C 

Ul 

o o o 

ooooooooo 



U. 

£ 


<x 

£ 


II 1 II 



£ 

d n 


a> 






Cl ffl 

© r 

CD 

_! 






o ▼ 

i r 

kO 

<1 






I o 

Ul < 

rs 

Z 






Ul -4 

<*• • 

d 

o 


d 0 '' 0 ©tr»rsrs©©* 4 ^© 



O' • 

•O d 



Z 

n O' 

▼ m — rsCD'orn-oo- 

d 


d © 

o» 



<r 

*» d < 

■O0 '«o©©rsd 0'0 


z 

d » 

m 

o 

M 

Ul 

m d 

-kOOO'O'O'O'COO' 

»- 

o 

▼ 

© 

i 

a 

£ 

• « * 

z 


«4 

• M 


z 


•4 H ^ 

«-*-4»-» 0000©0 


H- 

• H 

o *- 

« 

o 




o 

<z 

o »- 

z 


u 




£ 

£ 

1 z 

N Ul 

z 






►- 

■ UJ 

Ul 

Ul 


cn 



n 

z 

Ul « 

<_» o 

*4 


£ u 



-o 

Ul 

u a 

z - 

u 


o u 




u 

*n z« 

Ul £ 

mm 


z 



z 

£ 

rs u u. 

-1 u. 

U. 


£ <C O © 

0 

0 

0 

0 

0 

3 

9 

NriDONOKiriT'O 

3 

O 

-• -1 U. 3 UJ 

£ 


Ul £ 

♦ d CD 

'O^O'rsDdd^ 

OC 

u 

• 3 U. A O 

Ul 


A 3 

W 

•4 ^ 



O X O z u 

o 


£ O 



£ 

ID 

1 OC U 3 

u 


3 U 



O 

> 

• CD to cn 3 

t- z 

0 * 


z a 




o 

« 

H- 

<t 


OO'LIHZ O 
•OC»£C O U. 


k \ ■ _ 

<j * in a> 

o a 'Omh 

jzq^nz 

y <t Z * • li 

> <x 

£ 


0» , 

d cc 
* o 


£ O' 


O £ 
O 
£ 


£ in £ 

«£m 
lACZ 
:>> h 


< H K N ft. 

Jnnoo) 
ui -4 ro • z 
znno« 
AC • • I cc 

o © © *- 

O t 

CC -J 

<t u 

a ►» • ■ « 4 

oc ae oc oc to &. cc 

ZOOAZftd 

H u, u. uuuo 


<r r- . 
-I £ l 
li/ Ul ( 
OC £ ( 
QC O ( 
O £ < 


£ X 


£ 0 ©©©© 0 ©© 0 © 0 >. 0 » 0 * 0 ' 0 » 0 ' 0 '' 0 'C>' 0 ' 

OZOOOOOOOOOOO'O'O'O'O'OO'O'O'O- 
« 0 >CON'OIO»rtNHOO-M 1 fnoMD> 

<CMOOOOOOOOOOOOOOOOOOOO 

OC < I I I I I I I I » » 

K 3 

£>-111)11111 I I . I i | i l I 

u u 

UDOCDN'OrtfMN-OHf^MTrt'ONOO' 



OU.-4OOOOOOOOOOOOOOOOOOO 

u I I I I I I 


102 


END OP RUN 63 “POINT 

terminated; stop 


SAMPLE 62*22 

SUCCESSFUL (TEMP) RESTORE RUN43P22 AS (TP2 
CANCELLED: DDNANE RMDS UNKNOWN 
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END OF RUN 63 -POINT 

terminated; stop 


SAMPLE 63'26 

SUCCESSFUL (TEMP) RESTORE RUN63P26 

cancelled: ddname rmds unknqun 
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TERMINATED) STOP 


SAMPLE 63,2 7 

SUCCESSFUL < 1EMP ) RESTORE RUN63F27 AS <TP27> 
cancelled; ddnane rhds unknown 


o o N«naj’O>M-r)«a)'O 0 > 

0.** rKOONO-Kin^^TO 

w *- - 

Z <t ®HNOOOOO«---« 

« a » 


►- 

X 














X 


cs 



Cr- 

4 

o 

© 

CM 

« 

© 

o 

© 

© 

a 


4* 

•o 

oc. 

n 

m 

n 

m 

CM 

*« 

•9 

m 

© 

m 

Cl. 

u 

▼ 

o 

rs 

in 

M 

O' 

« 

rs 

K» 

M 

CM 

s 

rt 

Ul 



CQ 

4 

rs 


n 

r> 

V* 

si 

■o 

O 


o 

X 

u. 

n 

— 

© 

r> 


o 

o 

CM 

in 

in 

O 

© 

© 

« 

u. 


m 



• 









X 

U 1 

rj 


O 

o 

o 

© 

o 

o 

o 

o 

o 

o 

© 


o 

» 


1 

» 

1 


♦ 

1 

■ 

i 

« 

1 

» 


u 
























-•mo. 

4 © © 

a> cm 

n 

All 

© 

© 










4 CM *4 

cm © m 

n n 

« 

m 4 

rs 

o 









© 

ClrtT 

CM 4 CO 

m cm 

© 

© 4 

rs 

© 









X 

O' -J O' 

~ CO IS 

© © 

m 

m «« 

© 

© 









X 

o *• ** 

CM ** 

•4 mm 

-4 

^4 -4 

-4 

© 







(/» 

►- 

-1 

3 



© © © 

O O © 

© O 

© 

© © 

© 

© 







in 

Ul 

X 

Ui 


CM 4 4 

©40 

-• r» 

m 

© -4 

9 

IS 







o 

3 


nfcrt 

rs •* 9 

-* ^ 

o> 

o © 

4 

4 





© 


X 

•-» 

X 

co m 

CO 09 O 

m 4 

¥1 

© o 

O 

Cr. 





K> «* 



►- 

4 

msn 

N rt 4 

© rs 


s in 

c 

© 





© «* 


-1 

4 

Ui 

m C. ^ 

© o o 

© © 


o o 

© 

© 







X 

-1 

X 

• • • 

• * # 

♦ 4 

4 








0- X 


X 

— 


© © © 

0 

0 

0 

© © 

© 

© © 

o 

© 





X 


© 





1 

1 

* 


1 





CM — 

-a 















NOOH 

© 

-1 














o ▼ i n 

CM 

4 














• OP Ul © 

n 

X 














id “s rs, • 

CM 

O 



© C" <D 

CM © © 

4 n 

n 

© in 

o 

o 

N 




© - © CM 

4 

•-s 


X 

© © CM 

09 rs QQ 

4 n 

r» 

CM o 

O' 

© 

cm 




-4 © © 

4 

fr- 


4 

rs o r> 

> n n 

m 4 

o 

« © 

m 

m 


X 



*0 CM 

o 



Ui 

CM cm o 

O' O' CO 

o © 

m 

CD © 

9 

© 

*— 

o 



O' 4 

1 

© 


X 

♦ • 4 

4 • • 



• • 



X 

■4 



CM • ■ 


X 



•4 •* *4 

0 

0 

0 

© o 

o 

© © 

© 

© 


N- 



• ■ o *- 

• 

o 










o 

4 



o >- X 

►- 

3 










a. 

X 



X M Ul 

X 












►- 



l lil y ** 

Ul 



© 








n 

X 



U"UU 



U. 

Ul 








© 

Ul 



U U X *+ 

u 


o 

u 









u n 



X^ujIl 




Z 








X 

X CM 



yy jil 

u. 


a 

0 

0 

9 

0 

0 

0 

0 

V 

f\ © 4 

in m 

rs 

4 4 

n 

4000 

3 

a o 



jwau 

u. 


Id 

x 

4 O' 

cm in in 

O' O' 

© 

rs rs 

CM 

X 

u * 



3 ui m o 

Ui 


© 

3 


..—««« 







© 



©0X0 

o 


X 

o 








X 

m 



X O 3 

o 


3 

u 








O 

3 ■ 

<J cn cn 3 ►- X r*» 



X 

o 








U. 

© 

O X X ►- X o 4 

►- 












>- X 

O' 

X X 

O U. *• CM 

X 












K S 

■ 


U ~ O ►- CM 

o 











3 

*-» X 

<r 

« O' a *- < h h n 

a. 



— mm — 

^ ^4 «© ^ 

^ *4 ^ 

O' O' 

o 

o o 

O' 

O O' o o 

a. 

u 

x o» © 

<r ►- _i cs rs o 

© 


X 

© o o © o © o 

0©00'0'0'0'0' 

O' 

O' o O' o 


O o 


'•O'^Jiwn-o . 

X 


o 

xoooooooooo 

O' O' O' O' O' 

o 

O' O O' o 

3 

Z 

a o 

x id ui « n cm o 

4 



OOfflN 


cm n 4 

m 

© rs © O' 

a 

U 4 

X 


Ui X X X - * I 

X 


►- 









3 

4o©xxoooo 

*- 


4 

»-oooooooo©ooooooooooo 

« 

X 



0 0X0 • 



X 

4 111 

i i i i 

1 1 1 






»- 

X 

O 


X O X 

o 


u» 

3 








4 

-1 CM 

> 


XX 4 

u 


X 

►-III 

iiii 

1 • t 

1 1 

1 

* i 


I’ll 

O 

O' 

• t 

© a v- ►- • ■ © 

4 


Ui 

u 









4 Id 


xcaxxxxxcox 

X 


o 

u 

0 

9 

8 

7 

6 

3 

4 


m © rs 9 o 


M *« 

■ 

4 X 

~ 3 3 « X 3 

Ui 


X 










X ■ 

©axxxooaxx 

3 


o 

u. — o © 
1 1 1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


4 N 

t>> >*»-u.u,uuu 

O 


o 

1 1 1 f 

l i i 

1 






105 


END OF RUN 63 -POINT 27 
TERMINATED t STOP 


SAMPLE 63 1 28 

SUCCESSFUL < TEMP > RESTORE RUN63P28 AS (TF28) 
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END OF RUN 63 -POINT 28 

terminated: stop 


SAMPLE 63*29 

SUCCESSFUL (TEMP) RESTORE RUN63P29 AS <TP29> 
CANCELLED t DDNAflE RMDS UNKNOWN 
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END OF RUN 63 -POINT 29 

terminated; stop 


SAMPLE 63*30 

SUCCESSFUL (TEMP ) RESTORE RUN63P3Q AS (TP30> 
CANCELLED! DDNAMC RMDS UNKNOWN 
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END OF RUN 63 -POINT 30 

terminated: stop 


SAHPLt 63*31 

SUCCESSFUL fTCNP) RESTORE RUN63P31 AS UF3t> 
CANCELLED ! DDNAhE RttDS UNKNOWN 
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SAMPLE 63.32 

SUCCESSFUL (TEMF) RESTORE RUN63P32 AS <TF32> 
IANCELLED1 ODNAME RhDS UNKNOWN 
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END OF RUN 63 -FOINT 32 

terminated: stop 


SAMPLE 63*33 

SUCCESSFUL t TEMP ) RESTORE RUN 63 P 33 AS (TP 33 > 

cancelled: ddname rhds unknown 
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END OF RUN 63 -POINT 33 
TERMINATED J STOP 


SAMPLE 63*34 

SUCCESSFUL I TEMP ) RESTORE FUN63P34 AS (Tf3«> 
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END OF RUN 63 -POINT 34 

terminated: stop 


SAnPLL 31 »b 

SUCCESSFUL < TEflP ) RESTORE RUN 51 P 0 AS (TfQ) 
CANCELLED! I'DNANE RNDS UNKNOWN 
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L Nfi UF RUN 31 -POINT 

terminated: stop 


SArtFLL 'Hr? 

SUCCESSFUL (TEHM RESTORE RUNS 1 f 9 AS <TF9> 

cancelled: ddnane rhds unknown 
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END OF RUN 51 -P01N1 

terminated: stop 


SAMPLE SI ,10 

SUCCESSFUL f TEMP > RESTORE RUNS1M0 AS (TFlO) 

cancellid: ddnaiie kmos unknown 
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END OF RUN SI -POINT 10 

terminated; stot 


SAMPLE bit 11 
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terminated: stop 


snrtfLt ii.u 

SUCCESSFUL ( TEMP ) RESTOKE KUN51P12 AS < TF 1 
CANCELLED! DDNAME RMOS UNF NOUN 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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END OF RUN Si -POINT 
TERMINATED! STOP 
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END OF RUN 51 -POINT 14 

terminated: stop 


SAMPLE 51 » IS 

SUCCESS! Ul (TEMP/ RESTORE RUN 51 P 15 

cancelled: ddname rmos unknown 
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END OF RUN 51 -POINT 

terminated: stop 


SAMf LE til . 16 

SUCCESSFUL (TErtf) FESTORC RUN51P16 AS (TFI6) 
CANCELLED: DDNAME AMDS UNKNOUN 
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END OF RUN SI -POINT 16 

terminated; stop 


SAttfLC \*\,\7 

SUCCESSFUL <TEhP> RESTORE KUN51F17 AS < TF 17' 
CANCELLED: DDHA.lf HrtUS UNKNOUN 
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END OF RUN 31 -POINT 17 
TCMtlNATCM STOP 
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SUCCESSFUL ( TEttF ) ft ES.TORE RUN51F21 
CANCEL L f P 1 DDNAME KMDS UNKNOWN 
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END OF RUN 51 -FOINT 
TERMINATED 1 STOP 


SAMhLE 51* 27 

SUCCESSFUL (TErtf) RESTORE RUNMF27 
CANCELLED • DDNArtf RMDS UNKNCUN 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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END OF RUN 51 -POINT 

terminated; stop 
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* 


to *- 

z <x 

<£ tC 


O o> cd f i r« r ' v ^ ca f> - f i 
r m « ▼ f 1 oo < mi w « 

nciCD'r. 

®M-*OOOoOO^fi^ 


z 

«£ Ui 

o -* 
t u 
2: u. 
c u. 

JC t-u 
- Q 


n-o^fiO'Noono-f" 

ttNOOD(l«««(>"ONM» 

-TtHHOoommn 

rjooooooooooo 


X 

cc 


-ooO'CMfvroo-*<x.r^y)CK 

N'OnO^'OCDt>-"0<Dn 

"•OOOONH-f^Mn 

a)o>'On«Hf4O'®flD0e 

_'fM*4»4~4*4 V 4~4.«OO0 

000000000O00 


in 

c 


a 

o 


- z 

- <r 

« Ui 

-J c 
u 


NMMN'O-NCDlllta)** 
o^mOntcntino' *• o oo 

n.% 

-tCOC -00000000 

OCOOOOOOOOO 0 


Cl 

>- 

3 


o ® 

■ a: ^ 


VI 01 

a. Cl. 

c r 


O — I 


► ♦ o n 

> O o 


• • o * 

0 *- 

1 z ■ 1 
n ui uj * 
w - U ( 
our* 
Z*«uk 
Ui u. j 1 
JIlSI 
3 Ui CD < 
A O i£ I 
uC UD 

3 : 

1- X < 

a u. * 

U. *-• o » 


m Ac ▼ n o o ►* 


A O ▼ ►- . 


X O < 
<t X 


iX CD 

<c c 

» A OC 

: 3 z> 


U lC c 

> r cc o 
00c 
X u 

X 

o a k 
CC CC CC 
-• — 3 
X X o 
►“ *- u. 


ffl oc 
o o 
<t x n -h a. 
-/ OD Q3 O W 
Ui Pv ® . Z 

COK 4 O <t 
cc • • UL 
000 h- 
U N 

(X _J 

X <z o 

■ i A <t 
cc cc in cl ^ 
3 <r r 3 ui 
o a «x a. 3 
u. u o u o 


on'0O>9'^i'0n^'0CQ 
r» -«is>®^o»'0«-*rs^r 
MD-oNn^nOn^tfi 
~00000000<300 

OOOOOO0O0OOO 


OOO«K^ , 0»'O^^P'l®^fVMOOOOO 

TOM« NfflO^T 


to 

U. Ui 

o u 

X 

<t o 
— * CC 
au 3 
x u 
3 u 
X o 


-***•* — « , ^ ^ ^ ^ ^ O.CNO.O'O.O'fr'O'O-O- 
X 0000000000»»>CK(^^^»^0*9‘ 
3XOOOOOOOOOOO»C>*^0‘>OCKO'0^0 
h o > cd n ' 0 ^^nN-oo^(Nn»in>oNOO 



<£►-00000000000000000000 
CC « I I I I I I I • I I 

*- 3 

X t- I I t I I 1 1 1 t 1 . 1 1 1 1 1 I I 1 

ti o 

uaofflN'OkiTinMHO-N'nfn'OMDD' 
X 

au» — 0000000000000000000 


I 


I I 


I I I 


I 


124 


END OF RUN 51 -POINT 
terminated; stof 


SAMJ LE 31 . 29 

SUCCESSFUL l TEMP) RESTORE RUN51F29 

cancelled: ddname rmds unknown 
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END OF RUN 51 -POINT 

terminated: stop 


SAMPLE 5 1 » 30 

SUCCESSFUL ( TEMP' RESTORE RUN51F30 AS <TF30> 

cancelled: odname rmds unknown 
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END OF RUN 51 -FOINT 30 

terminated: stop 


SAMPLE 01,31 

SUCCESSFUL (TEMP) FES T ORE RUN51P31 AS (TF31) 
CANCELLED < DDNAME RMD5 UNKNOWN 
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END OF RUN SI -POINT 31 

terminated; siop 


SAMPLE SI. 32 

SUCCESSFUL (TEMP) RESTORE RUN51P32 AS <TF3 
CANCELLED: DDNAME RMDS UNKNOWN 
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END OF RUN 31 -POINT 
TERMINATED l STOP 


RADIAL VELOCITY VS CONCENTRATION 
Z-152 NM AND R/R0--0.275 
NO- 999 AND N«-990 
VDAR- 0.0343 MPS 

VRMS- 0. 177A NPS 


ORIGINAL PAGE IS 
QF POOR 0 M 4 LfTY 


o o 

Cl. — 
(Si ►* 

x <Z 
<c cc 


— non »-on®noOP> 

rt ® O N rt f l M*< w ♦ 
OO OO 


<C u. 
O -I 

a. w 

•r. ** 

Z u. 
<t u. 
CC UJ 

w o 
u 


--on-flO-oo^o^o 
rv^mrv»«rfs.(K-o>o)0^ 
rv^3ors.r^o(Oa3C^-^-*© 
« *• r> w ~ <* <► m in *■ >o 

NOON>N'0^^t'0<Q 
OOOOO - NfO T-Ortlrt 


rvrt^c<w*r>oooor‘0 

'OtO.N«Ot’rinO(CO 

oooooooooooo 


m (rt 
c a 4. 

y. r 

s 

r* io 
r> w o o 

0 is. i oc 

1 » y •< 

y o n * 
o * o w 
-* O «0 
O I - 
▼ M 

n - ■ 
- « o — 


cK>^<«o-rsir, P ‘>.rv<rv 
»nf •.no-ONO^rt'O 

m(N(Nin-cMO-®oo^ffl 

OOOO — 

oooooooooooo 


N>fflWN^O<«ONO 
► MOON-OOOHrtOO. 
HOKl».TnrtM»NN« 

oooo**-Wf<r<rt-'<s 

OOOOOOOOOOOO 


y M uu 
u <J * •* 
x — ui u. 
y U. j u. 
JU.3U 
3 UJ A O 

a o a u 

OC U 3 


u. — a 

O ' 


o r* *— 

•- n oc 

*- NO 

< n n m a. 


c h j > o o m 

►- j x ui r» © • x 

ZUUQC «<«<04 
UJ oc X X • • cc 

X oc a Q o o ►— 

o o c U ■ 

XU oc -l 

XX < ^ 

q a ■ ■ a <z 
QCKCC«0C<A(LC 
«*-* 33 <IX 3 uj 
XXOOOXO.O 
►“►“U.U.UOUO 


© 

U. UJ 
O U 

z 

^«00000000"»«n04NM()N) 

u *< HMrto.nr»NH 


3 U 

X o 


Z 0o000000000*9»0'>0»0>0*9>9‘0o 
c ZOOOOOOOOOOO-O'O-OO'OOOOCK 

«o>fflNon»np<«oo*'(NrttH’ON®> 

<*-00000000000000000000 

X « I I < I t I I I • I 

*- 3 

X *- » I • I I I I i I I I I ■ I I I I I I I 

UJ o 

U 3 O O* O N TiONnO^Nfn vn -ON® > 

Z J 

Ou. — ooooooooooooooooooo 

U • I I I I t I | I I | 


129 


END OF RUN 31 -POINT 33 

terminated; stop 


SAMPLE 51, J4 

SUCCESSFUL (TEMP) KESTORE RUN51F34 AS iTF34> 
CANCELLED. DDNAHE RMDS UNKNOWN 
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END OF RUN 51 -POINT 34 

terminated: stop 


SAMPLE 31,35 

SUCCESSFUL < TENF ) RESTORE RUN51F35 

cancelled: hdnami rnds unknown 
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